Pure spin current in a two-dimensional topological insulator 
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We predict a mechanism to generate a pure spin current in a two-dimensional topological insulator. 
As the magnetic impurities exist on one of edges of the two-dimensional topological insulator, a gap 
is opened in the corresponding gapless edge states but another pair of gapless edge states with 
opposite spin are still protected by the time-reversal symmetry. So the conductance plateaus with 
the half-integer values e 2 /h can be obtained in the gap induced by magnetic impurities, which means 
that the pure spin current can be induced in the sample. We also hnd that the pure spin current 
is insensitive to weak disorder. The mechanism to generate pure spin currents is generalized for 
two-dimensional topological insulators. 
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Since graphene, a single-layer honeycomb lattice of 
carbon atoms, has been prepared laboratorially by 
Novoselov et al.,[l| it has attracted considerable atten- 
tions due to its novel properties in condensed matter 
physics and potential applications in devices. [2H10I] 
Graphene is the first independent two-dimensional (2D) 
crystal that has been experimentally achieved, which 
leads to new interests in 2D systems. For example, the 
Kane-Mele model, a quantum spin Hall effect (QSHE) 
was first proposed in graphene with spin-orbital coupling, 
which is the first example of topological insulator. [Il|,|l2j 
Topological insulators are time-reversal symmetric sys- 
tems whose intrinsic spin-orbit coupling (SOC) opens a 
bulk gap while generating the Kramers doublet of edge 
states owing to the nontrivial Z2 invariants of the occu- 
pied bands. The edge states force electrons with oppo- 
site spin to flow in opposite directions along the edges 
of the sample, which lead to quantized spin Hall conduc- 
tance. However, the intrinsic SOC in realistic graphene 
is quite weak and the gap opening was small, so the 
QSHE in graphene is difficult to be observed. [13J Never- 
theless, recently a monolayer honeycomb lattice of silicon 
called silicene has been synthesized and attracts much 
attention. [l44l7| Silicene has a relatively large intrinsic 
spin-orbit gap of 1.55meV, as makes experimentally ac- 
cessible the Kanc-Mclc type QSHE. 16] 

The development of the topological insulator opens a 
new and powerful way for the spintronic applications 
due to its spin-dependent edge states. How to gener- 
ate pure spin currents in low-dimensional systems is the 
main challenge in the field of spintronics. The aim of 
this work is to propose a method for generating a pure 
spin current in a 2DTI. In this Letter, as a concrete ex- 
ample, we theoretically study the electron transport in 
Kane-Mele model with magnetic doping at one edge, as 
shown in Fig. 1. Most of the results are also applicable 
to general 2DTIs. Before presenting our detailed calcula- 
tions, we first analyze why pure spin current can be gen- 



erated in the present device. The intrinsic SOC which 
originates from intra-atomic SOC, converts the sample 
into a topological insulator with a QSHE.[ll| The gap- 
less edge states are protected by time-reversal symmetry 
and is thus robust to non-magnetic impurities that do 
not break this symmetry. But a pair of edge states are 
destroyed when magnetic impurities exist on this corre- 
sponding edge (see Fig. 1). Because another pair of edge 
states, with opposite spins containing opposite propaga- 
tion directions, are still protected by time-reversal sym- 
metry, we can observe a pure spin current in the sample, 
which is confirmed by our following calculations. 

In the tight-binding representation, we consider 
the Kane-Mele Hamiltonian defined on a honeycomb 
lattice: [H [l2 1 



H = t J2 c L c 3<?+ iX J2 V V C L S z C 3<r 

+ ia Yl C L( S x d ^)zc 3 a + K^2^c\ a c lrT . (1) 

The symbols (ij) and ((ij)) denote the nearest and the 
next-nearest neighbors, respectively, and a ='t ) l (or ±1) 
denotes spin index. The first term is the nearest-neighbor 
hopping. The second term describes the intrinsic SOC. 

va A is 



Here the site-dependent Haldane phase factor [11 1 
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defined as i/y = (di x d2)/|di x d.2 1 = ±1, where dj 
denotes the vector from one atom to one of its nearest 
neighbors. s z is a Pauli matrix describing the electron's 
spin. The third term is a nearest neighbor Rashba SOC 
term, which can be produced by applying an electric field 
perpendicular to the sheet. The fourth term is a stag- 
gered sublattice potential (& = ±1). It is interesting to 
notice that, Eq. (1) is almost applicable to silicene ex- 
cept for the Rashba SOC term, which is present between 
the next-nearest neighbors in silicene but between the 
nearest neighbors in our model. The focus of this work is 
the introduction of magnetic impurities (red square dots 



in Fig. 1) on the uppermost zigzag chain of the sample 
(n = 1), 



0.5 






(2) 



where M is the strength of exchange interaction induced 
by the magnetic impurities. This term breaks the local 
time reversal symmetry on the upper edge. 




FIG. 1: (Color online) Schematic diagram of a zigzag honey- 
comb lattice nanoribbon with magnetic impurities (red square 
dots) on the upper edge. Its unit cell is marked by the dashed 
lines. The red and blue arrows at the edges denote propaga- 
tion directions of the opposite spins in the edge states. 

In the following numerical calculations, we use the hop- 
ping energy t as the energy unit. The width N is chosen 
as N = 50 in all calculations and the nearest neighbor 
atom-atom distance is a. The strengths of the intrin- 
sic SOC, the Rashba SOC and the staggered sublattice 
potential are A = 0.06£, a — 0.05t, and X u — O.lt, re- 
spectively. These parameters define the system as a two- 
dimensional topological insulator, [llj 

First, we investigate the energy subbands obtained by 
solving the lattice model in strip geometry, as shown 
in Fig. 2. In the pure case (M = 0, Fig. 2(a)), the 
edge states traverse the energy gap in pairs. The gap- 
less edge states are robust against small non-magnetic 
perturbations since they are protected by time reversal 
symmetry. [Ill] However, in the presence of magnetic im- 
purities (M ^ 0, Fig. 2(b)-(d)) at the upper edge, the 
corresponding pair of gapless edge states is destroyed 
and a gap can be opened due to the local time rever- 
sal symmetry breaking. Moreover, the magnitude of the 
gap opened by the magnetic impurities increases with en- 
hanced M, and can reach O.lSey for M = l.Ot. Another 
pair of gapless edge states, still protected by time reversal 
symmetry, persist on that edge without magnetic impu- 
rities. But the electron-hole symmetry is broken in the 
preserved gapless edge states, which cross at ka — 7r even 
for a very small value of M . 

Next, we examine the influence of the magnetic 
impurities on the conductance of the system, as 
shown in Fig. 3. The two-terminal conductance of 
the system can be calculated by the nonequilibrium 
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FIG. 2: (Color online) Calculated energy bands in the zigzag 
honeycomb lattice nanoribbon in the presence of the intrinsic 
spin-orbit coupling and the Rashba spin-orbit coupling for 
M = (a), M = 0.4t (b), M = 0.8i (c), and M = l.Ot (d). 
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3: (color online) The conductance G vs E for different 



Green's function method and the Landauer-Biittiker for- 
mula G{E) = ^Tv[T L (E)G r (E)T R (E)G a (E)}, where 
T P (E) = i\Pp{E) - ££(£)] is the line- width function 
and G r (E) = [G a (£)]t = 1/[E - H ce „ - E£ - ££] is 
the retarded Green function with the Hamiltonian in the 
center region H ce n.[18| The self-energy ££ due to the 
semi- infinite lead-p can be calculated numerically. [19( In 
the case of M = 0, the quantized conductance plateau 
appears, with the plateau value 2e 2 /h coming from the 
contributions of two pairs of the gapless edge states. For 
greater values of M, the conductance plateau 2e 2 /h is 
suppressed and evolves into a conductance plateau e 2 /h 
in the energy gap opened by the magnetic impurities. 
The conductance plateau e 2 /h is only contributed from 
the gapless edge states at the lower edge without mag- 



netic impurities, which can be seen from Fig. 4. Without 
magnetic impurities (M = 0), there are two pairs of per- 
fect edge states with opposite spins at two edges of the 
sample, so the currents through the sample is spin unpo- 
larized. However, when the gapless edge states at the up- 
per edge are destroyed by magnetic impurities (M ^ 0), 
the spin- velocity locked channels persist only at the lower 
edge of the sample. Therefore, the current in the sample 
consists of the opposite spins moving in opposite direc- 
tions that is called pure spin current, which is just the 
aim of the present work. For a definite arrangement of 
bias voltage, there remains only, say, right going channels 
with spin up working. Moreover, the pure spin current 
in the sample is really invariant under time reversal. 



The spin-resolved conductance matrix can be written as 






(4) 



which can also be calculated by generalized Landauer 
formula for spin transport. The conductance Gff and 
Gfi can be obtained when we assume that only spin-up 
electrons are injected from the left lead into the sample 
and collected in the right lead. We can also calculate 
Gif and Gff in the same way by assuming only spin- 
down electrons are injected from the left lead. The total 
conductance G and the spin polarization P in lead-i? can 
be respectively defined as 20, 21] 



G — G^ + Gff + Gff + G 



-14- 



(5) 



and 



to 

o 
3 o 



o 

Q 



1 


(a) 






J 




left 

right 










// ' 





100 




100 



FIG. 4: (Color online) Spin-resolved local density of states 
(LDOS) along the width of sample for (a) M — and (b) 
M = 1.0. The energy is chosen E = 0.02t in the simulation. 
N a — 2n — 1 and N a = 2n correspond to spin-up and spin- 
down LDOS, respectively, where N a - 1,2, •■• ,2N. The 
black up and down arrows denote spin up and spin down. 
The blue and red curves denote left and right propagation 
direction 
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To test the above arguments in a more direct way, we 
also studied the spin-resolved conductance and spin po- 
larization when the magnetic impurities exist on the up- 
per edge of the sample. For the sake of simplicity, we 
assume that in the leads, the staggered sublattice poten- 
tial, the intrinsic and Rashba SOC do not exist, i.e., the 
Hamiltonian of lead- p is simply 






(3) 



FIG. 5: (Color online) Spin-resolved conductance (a) and spin 
polarization (b) vs E for M — 1.0. The parameters A = 0, 
a — 0, and \ v — are chosen in leads 

Fig. 5 (a) and (b) show the spin-resolved conductance 
G and spin polarization P versus the energy E for M = 1. 
In Fig. 5 (a), the total conductance manifests itself with 
the plateau value e 2 /h due to the presence of the mag- 
netic impurities on the upper edge of the central region. 
Due to the topological nature of the edge states, this 



plateau is insensitive to the mismatch between the sam- 
ple and the leads. We can also find that the spin-up and 
spin-down electrons are not mixed when they transport 
through the sample, i.e., G-f-j, = G^ = 0. The spin polar- 
ization can almost reach 100% in the energy gap opened 
by the magnetic impurities [see Fig. 5 (b)] because the 
spin-up electron can fully transport through the sam- 
ple while the spin-down electron can hardly transport 
through the sample in the gap [see Fig. 5 (a)]. Beyond 
the gap, due to the conduction band mismatch between 
the central region and the leads, there are resonant tun- 
neling peaks in the conductance beyond the gap induced 
by the magnetic impurities, as shown in Fig. 3 and Fig. 
5. 




FIG. 6: The conductance G as a function of the energy E 
for the disorder strengths w = 0.5£ (a) and as a function of 
the disorder strengths w for the energy E = 0.02t (b). The 
error bars show standard deviation of the conductance for 100 
samples. 

Finally, we examine the non-magnetic disorder effect 
on this spin-polarized conductance plateau e 2 /h. Due 
to disorder, random on-site potential Wi is added for 
each site i in the central region, where Wi is uniformly 
distributed in the range [—w/2,w/2] with the disorder 
strength w. Fig. 6 (a) and (b) show the conductance G 
versus the energy E at the disorder strength w = 0.5t and 
G versus the disorder strength w at the energy E — 0.02£, 
respectively. The results show that the quantum plateau 
of e 2 jh is very robust against non-magnetic disorder be- 
cause of the topological origin of the edge states. The 



quantum plateau maintains its quantized value very well 
even when w reaches l.Ot. The robust and stable plateau 
of e 2 /h means that the pure spin current of the system 
is insensitive to weak disorder and protected by time- 
reversal symmetry. In addition, even for a large disor- 
der strength w (e.g., from w = l.Ot to w = 3.0t), the 
conductance is increased rather than decreasing with the 
increasing disorder strength. This is because although 
the strong disorder weakens the edge states, it also result 
in the mobility of the energy band structure, 22j so the 
value of G increases in the range of w = l.Ot to w = 3. Of. 
With further increasing of the disorder strength, the con- 
ductance gradually reduce to zero, the system eventually 
enters the insulating regime. 

In summary, we predict a new mechanism to generate 
a pure spin current in a two-dimensional topological in- 
sulator. As the magnetic impurities exist on one edge of 
the sample, the corresponding gapless edge states is de- 
stroyed but another pair of gapless edge states with op- 
posite spin are protected by time-reversal symmetry. So 
a pure spin current with the spin-up and spin-down carri- 
ers moving in opposite directions can be observed in the 
system. Moreover, the pure spin current has also been 
found to be robust against non-magnetic disorder. The 
mechanism to generate pure spin currents can be gener- 
alized for two-dimensional topological insulators, such as 
HgTe/CdTe quantum well and silicene nanoribbons. 
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